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Abstract 
(N' 

We obtain the binding energy of an infinitely heavy quark-antiquark pair from 
Dirac brackets by computing the expectation value of the pure QCD Hamilto- 
nian. This procedure exploits the rich structure of the dressing around static 

o. 

CO , fermions. Some subtle points related to exhibing explicitly the interquark 

O 

t> 

o 



energy are considered. 
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It is well known that to gain a deeper insight into gauge theories a proper study of the 

concepts of screening and confinement is of considerable importance. In this connection 

the binding energy of an infinitely heavy quark-antiquark pair represents a fundamental 

concept in QCD which is expected to play an important role in the understanding of quark 

confinement. In this respect we recall that asymptotic freedom allows one to use perturbation 

theory for the description of high energy phenomena. But within this framework, we cannot 

explain low energy phenomena such as the permanent confinement of quarks and gluons. The 

reason is that the infrared divergences and gauge dependence make bound-state equations 

very hard to approximate. We further note that the choice of the gauge has a strong 
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influence on the propagator. For instance, the studies of Caracciolo et al. H] were crucial in 
the investigation of the ambiguity in the definition of the gluon propagator. 

In view of the above-mentioned difficulties much attention has been recently devoted 
to formulations of QCD in which gauge-invariant variables are explicitly constructed. By 
doing so it has been possible to obtain a more direct physical insight into the description 
of charged fields P|§. As a resu h of this development the matter fields (quarks) are now 
dressed by a cloud of gauge fields. Using this approach, the pure QCD correction to the 
Coulomb potential in 3+1 dimensions has also been derived more recently by Lavelle et al. 
||. In particular, it was showed that such correction contains a dominant anti-screening 
contribution and another one which corresponds to screening by physical gluons. 

On the other hand, it is well known that gauge theories fall into the class of constrained 
systems, where a systematic procedure for quantizing such systems has been given by Dirac 
0. The point we wish to emphasize, however, is that because of the structure of the 
constraints the resulting equal-time commutation relations involve in general the coupling 
constant. This raises the question of how to recover from them the interquark potential 
between charged fields following the conventional path via the expectation value of the 
QCD Hamiltonian. This problem is addressed in the present letter. 

II. INTERQUARK POTENTIAL 

Before computing explicitly the interquark potential, we shall first present the Hamilto- 
nian analysis for the Yang-Mills field coupled to an external source J°. We start with the 
Lagrangian 

£ = --Tr (F^F^) - A a J° = —F^F 4 ^ - A a J°. (1) 

Here A^ (x) = A a (x) T a , where T a is a hermitian representation of the semi-simple and 
compact gauge group; and F" = d^A^ — d u A a ^+ g f abc A b ^A c u , with f abc the structure constants 
of the gauge group. The Dirac procedure |5| as applied to (|I|) is straightforward. The 



canonical momenta are U afl = —F a0fJ- , which results in the usual primary constraints Ilg = 0, 
and U ai = F ai °. The canonical Hamiltonian following from the above Lagrangian is: 

H c = f d 3 x (--W}U ai + HftfA 00 + jF^F aij - gf abc Ii ai A m A1 + A a Q J°) . (2) 

The persistence of the primary constraints leads to the following secondary constraints 

Q a W ( x ) = diU m + gf abc A bi Ii1 - J° « 0. (3) 

It is easy to check that there are no further constraints in the theory, and that the above 
constraints are first class. The corresponding total (first class) Hamiltonian that generates 
the time evolution of the dynamical variables is given by 

H = H c + J d 3 x (eg 0) U a (x) + c\ 0) fi a(1) (x) ) , (4) 

where Cq and c" are arbitrary functions. Since IIq p» for all time and Aq (x) = [A% (x) , H] = 
Cq (x) , which is completely arbitrary, we discard Aq (x) and IIq (x) because they add nothing 
to the description of the system. The Hamiltonian then takes the form 

H= f d 3 x (--U^U ai + jF^F aij + c a (x) [dmi + gf^A^Wr -J )), (5) 

where c a (x) = c\ (x) — Aq (x). 

Therefore, we have the first class constraints fW 1 ) (x) , which appear at the secondary 
level. Now the presence of the arbitrary quantities c a (x) are undesirable since we have no way 
of giving them a meaning in a quantum theory. To circumvent this trouble, we introduce a 
supplementary condition on the vector potential such that the full set of constraints becomes 
second class. A particularly appealing and useful choice is given by 
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Jd\{x-i) k At ] {i + \{x-i))^^ (6) 

b 



where A(0 < A < 1) is the parameter describing the spacelike straight path x k = £ fc + X(x — 
O', on a fixed time slice. This supplementary condition is the non-Abelian generalization of 
the gauge condition discussed in |||7|], which leads to the Poincare gauge ||. For simplicity 
we restrict our considerations to £ fc = 0. As a consequence, (§) becomes 



1 



fiP) (x) = / d\x k A a k (Xx) « 0. (7) 

o 

Now we come to the calculation of the Dirac brackets. By following the Dirac's procedure 
one arives at 

{AUx),Ai(y)Y = = {uUx) 1 Ili(y)}\ (8) 

{At (x) , Tl b > (y)}* = 8 ab S{S^ (x - y) - J dX U ab -^- - gf ahc A\ (x)) x^^ (Xx - y) . (9) 

Note the presence of the last term on the right-hand side which depends on g. In passing 
we note that similar Dirac brackets were obtained independently in reference [0]. 

We are now equipped to compute the interaction energy between pointlike sources in 
pure QCD, where a fermion is localized at the origin and an antifermion at y. In order 
to accomplish this purpose, we will calculate the expectation value of the energy operator 
H in the physical state \Q), which we will denote by (H)q . From our above discussion, we 
see that {H)q reads 

(H) a = (n\Jd 5 x (^n«n» + \f%f*) |n> . (10) 

Since the fermions are taken to be infinitely massive (static), this can be further simplified 



as 



(H) n = (o| J d'x (-In? ( X ) n« (*)) |n> . (n) 

Let us also mention here that, as was first established by Dirac [|H]], the physical states 
\Q) correspond to the gauge invariant ones. It is helpful to recall at this stage that in the 
Abelian case \Q) may be written as |J 

\Q) = |*(y)*(0))=^(y)exp Ugjdz^A, (z) J i> (0) |0) , (12) 

where |0) is the physical vacuum state and the line integral appearing in the above expression 
is along a spacelike path starting at and ending at y, on a fixed time slice. It should be 



clear from this discussion that the strings between fermions have been introduced in order 
to have a gauge-invariant function \Q). According to this viewpoint the fermion fields are 
now dressed by a cloud of gauge fields. 

The expression (|T2| ) may be extended on account of the fact that we have non-Abelian 
fields. Accordingly, we can write a state which has a fermion at and an antifermion at y 



as 



|fi)=^(y)C/(y,O)V(O)|0>, (13) 

where 

U(y, 0) = Pexplig I dz l A* (z) T a J . (14) 

As before, the line integral is along a spacelike path on a fixed time slice, P is the path- 
ordering prescription and |0) is the physical vacuum state. 

This last point gives us an opportunity to compare our work with the standard Wilson 
loop procedure |11|| , to make sure that the known results are recovered from the general 
expression flilf) in the weak coupling limit. In effect, due to asymptotic freedom, the short 
distance behavior of the interquark potential is determined by perturbation theory. Accord- 
ing to this, at weak coupling, one can expand 

*7(y,0) = Pexp lig j dz* A* {z) T a \ = P (l + igfdz i A^(z)T a + ...\ . (15) 

This implies that, at lowest order in g, the non-Abelian generalization of the dressing frame- 
work is the same as in the Abelian theory. Thus, at short distances, one should expect to 
obtain the known Coulomb potential in addition to a correction of order g 4 , as we will now 
show. 

With this in view, our next task is the computation of the expectation value of H in the 
physical state \fl) given by the expression (pTT|) . From the above Hamiltonian analysis one 
distinguishes here an Abelian part (proportional to Cf) and a non-Abelian part (propor- 
tional to the combination CfCa)- We first consider the Abelian part which is identical to 
the QED case. To do this, we shall begin by observing that 



y 
II? (x) \Q) = $ (y) U(y, 0)V (0) II? (x) |0) + ^T a J d Zl 5 (x - z) \Q) . (16) 

b 

Using this in (|TT|) we then evaluate the interaction energy in the presence of the static 

charges 

y y 
V9 2 = - g HrT a T a J dz l f dz[5 (z - z') , (17) 



remembering that the integrals over z l and z[ are zero except on the contour of integration. 
Writing the purely group theoretic factor trT a T a = Cf, the expression ( |TTD leads to 

V°\L) = l -g 2 C F kL, (18) 

after subtracting the term (H) = (0\H\0), where |y| = L and k = 5^(0). This calculation 
shows that special care is required in order to clarify the appearance of this peculiar result. 
It may be remarked, however, that the origin of the divergence is quite clear, so that it 
is possible to extract the Coulomb potential from the infinite contribution. Notice that 
the origin of the divergent factor k is due to the fact that the thickness of the string is 
nonvanishing only on the contour of integration. It is worth stressing at this stage that a 
more careful examination of the term ^- / d 3 x J dziS^ (x — z) ] reproduces exactly the 
expected Coulomb interaction between charges after subtracting the self-energy term, as was 
discussed in ||. Having made this observation, we write immediately the standard result 
for the potential to order g 2 , that is, 

Let us also mention here that if we had considered a modified form for the supplementary 
condition (§), which is equivalent to the Coulomb gauge 0, the result for the potential 
would have been the same. 

We now turn our attention to the non-Abelian part. From the expressions @ and (|TT1), 
we see that the g 4 contribution may be written in the form 



v 9 ' = d 3 x(o\(r) |o), (20) 



where 

y i 

I* = g 2fac T b I dz k I ^^ ^ z i 5 ^ _ Az ^ _ ^ 21 ^ 



The expression (^T[) may then be further manipulated as described in [[]]. Thus, if we use 

spherical coordinates, we find that 

y 



p = g 2 fb ac Tb * J dz k A c {z) £ y* y, ^ y^ ^ ) (22) 

z xr •/ i 



Now, by employing (^) we can reduce (p0|) to 

^ 4 = ^r / d 3 x <0| ( /> ac T b f--±- 2 1 dz k At(z) ^ YL (0', <<*) Y lm (6, <p)) |0) , (23) 

\ I I I I Ira ) 

which, by introducing the integration variable r = x and using usual properties for the 

spherical harmonics, may be rewritten as 

y y 
V 9 " (L) = -q^CaCf (-jjjdz'Jdz'W^z'). (24) 



Here Dij(z,z f ) is the propagator, which is diagonal in colour and taken in an arbitrary 
gauge. Thus, in order to carry out this calculation, we choose, for example, Dtj(z, z') in the 
Feynman gauge. Hence expression ( Pip reduces to 

V* (L) = -J^CxC, (_ 1) /dz/dz'^. (25) 

o o v ; 

This allows us to derive the g 4 contribution 

V* 4 (L) = -g*±- 2 C A C F ^ log {nL) , (26) 



where \i is a cutoff. By putting together Eqs.(|19D and (|26|), we obtain for the total interquark 
potential 

y (L) = -9'Cf^ (l + ^C A log (fiL)) . (27) 

In this way one obtains the known heavy interquark potential at lowest order in g [Q. 
However, the central difference between the above analysis and that of Ref. |4j] rests in the 



fact that the potential (p7| ) is directly recovered from the constraints structure of the theory 
we have discussed. In this context, the present investigation complements the discussion 
done in Ref. ||, as well as it reveals the general viability of our analysis. Thus it seems 
a challenging job to extend the scope of applicability of the above analysis. We expect to 
report on progress along these lines soon. 
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